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We consider Ising models defined on periodic approximants of aperiodic graphs. 
The model contains only a single coupling constant and no magnetic field, so the 
aperiodicity is entirely given by the different local environments of neighbours in 
the aperiodic graph. In this case, the partition function zeros in the temperature 
variable, also known as the Fisher zeros, can be calculated by diagonalisation 
of finite matrices. We present the partition function zero patterns for periodic 
approximants of the Penrose and the Ammann-Beenker tiling, and derive precise 
estimates of the critical temperatures. 



1 Introduction 

The Ising model is the paradigm for a second-order phase transition in two- 
dimensional models of statistical mechanics. Without an external magnetic 
field, the partition function of the two-dimensional Ising model can be calcu- 
lated explicitly for regular lattices. However, if the Ising model is defined on 
an aperiodic graph, for instance on the Penrose tiling (PT), this is no longer 
the case. Therefore, the influence of an aperiodic order in the underlying 
structure on the phase transition of the Ising model has been investigated by 
various means, including numerical simulations, series expansions, and zeros 
of the partition function; an overview on the results and a comprehensive lists 
of references on the subject can be found in recent review articles. 11113 

The most general prediction stems from heuristic scaling arguments, 
adapted from a relevance criterion for disordered Ising modelsB It yields an 
inequality involving a characteristic exponent that describes the fluctuations 
in the disordered or aperiodically ordered system, the correlation critical ex- 
ponent v of the pure system, and the space dimension of the fluctuation. 
Planar quasiperiodic graphs obtained by cut-and-project methodsu have low 
fluctuations, because they are flat sections through higher-dimensional peri- 



Figure 1. The first three periodic approximants (m = 1, 2, 3) of the Ammann-Beenker tiling. 



odic lattices. Therefore, the aperiodicity is expected to be irrelevant in these 
cases, which conforms with the results obtained for specific examples, notably 
the PT. 

Here, we consider partition function zeros obtained for periodic approxi- 
mants of the rhombic PT and the Ammann-Beenker tiling (ABT) . The unit 
cells of the first three approximants for the latter are shown in Fig. [j]. 

2 The Ising model 

We place Ising spins a- k e { — 1,1} on the vertices of the periodic approxi- 
mants. Spins interact by a nearest-neighbour interaction J. The energy of a 
configuration a — {a ■ ■ | 1 < j < M} G { — 1,1} on a finite graph with M 
vertices is 

where the summation is performed over all pairs (j, k) of neighbouring spins, 
i.e., those located on vertices that are connected by an edge. The correspond- 
ing partition function is the sum on all configurations 

Z(/3) = ]Texp[-/?£(cT)], (2) 

a 

where (3 is the inverse temperature. 
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Figure 2. Partition function zeros of the square lattice Ising model. 



3 Partition function zeros 

We are interested in the pattern of zeros of this function in the complex tem- 
perature plane. In the thermodynamic limit M — > oo, the zeros accumulate, 
filling curves or areas in the complex plane. These separate different regions 
of analyticity of the free energy, which is essentially the logarithm of the 
partition function, and hence correspond to phase transitions. 

The zero pattern for the square lattice case is shown in Fig. ^, which 
displays the zeros in the complex variable z = exp(2/3J). It consists of two 
circles. The two intersections at z — \/2 + 1 and z = \/2 — 1 with the positive 
real axis correspond to the ferromagnetic and antiferromagnetic critical points 
w c = tanh(/5 c J) = ±\/2, respectively. 

4 Results and conclusions 

The partition function zeros of the Ising model on any periodic planar graph 
with N spins in a unit cell caEL.be calculated by diagonalising a two-parameter 
family of AN x AN matrices.Qcl The powerful tool behind this method is the 
Kac-Ward determinant or Pfaffian method fbr calculating the Ising model 
partition function on general planar graphs.13 Thus, this method allows to 
compute partition function zeros for Ising models on infinite periodic graphs 
with an arbitrary numerical precision. 

The result for the first approximants of the PT and the ABT are shown 
in Figs. ^| and ^. These show accumulated zeros obtained by diagonalising 
the Kac-Ward matrix for various values of the two parameters, again in the 
complex variable z = exp(2/3J). The corresponding values of the critical 
temperature are listed in Table These yield precise estimates for the critical 
temperature for the aperiodic tilings, the errors were estimated on the basis 
of the apparent convergence of the first few terms. 
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Figure 3. Part of partition function zeros for periodic approximants of the Penrose tiling. 



Figure 4. Partition function zeros for periodic approximants of the Ammann-Bccnkcr tiling. 



The zero patterns, besides the trivial symmetry under complex conju- 
gation of z, still show the z <-» z^ 1 symmetry that is a consequence of the 
bipartiteness of the tilings. This also implies that the ferromagnetic and an- 
tifcrromagnctic phase transitions are of the same type. However, some of the 
symmetry of the square lattice case, corresponding to the self-duality of the 
square-lattice Ising model, is lost. In fact, the dual graphs of the aperiodic 
graphs under consideration contain triangles, hence are not bipartite. The 
antifcrromagnetic model on these dual graphs will suffer frustration. This 
may be responsible for the rather complicated, presumably fractal or frac- 
tally bounded zero patterns, particularly further away from the positive real 



Table 1. Critical temperatures w c = tanh(/3 c J) for periodic approximants of the PT and 
the ABT. Here, m labels the approximants with N spins per unit cell. 



m 


Penrose tiling 

N w c 


Ammann-Beenker tiling 

N w c 


1 


11 


0.401440 380 


7 


0.396 850 570 


2 


29 


0.395 411099 


41 


0.396 003 524 


3 


76 


0.395 082 894 


239 


0.395 985 346 


4 


199 


0.394 554 945 


1393 


0.395 984 811 


5 


521 


0.394 523 576 


8119 


0.395 984 795 


6 


1364 


0.394454 880 


47321 


0.395 984 795 


7 


3571 


0.394451035 






8 


9349 


0.394441450 






9 


24476 


0.394439 826 






10 


64079 


0.394439 319 






oo 


oo 


0.394 439(1) 


oo 


0.395 984 79(1) 



axis. The zeros around the two phase transition points still lie on a circle 
intersecting the real axis orthogonally, consistent with the fact that the phase 
transition belongs to the same universality class as on the square lattice. 
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